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ABSTRACT 

We present deep 3.6-8/im imaging of the Hubble Deep Field South with IRAC on the Spitzer Space Tele- 
scope. We study Distant Red Galaxies (DRGs) at z > 2 selected by J s -K s > 2.3 and compare them to a sample 
of Lyman Break Galaxies (LBGs) at z = 2— 3. The observed UV-to-8/im spectral energy distributions are fit 
with stellar population models to constrain star formation histories and derive stellar masses. We find that 
70% of the DRGs are best described by dust-reddened star forming models and 30% are very well fit with old 
and "dead" models. Using only the I— K s and K s -4.5pm colors we can effectively separate the two groups. 
The dead systems are among the most massive at z ~ 2.5 (mean stellar mass < M* >= 0.8 x 1O U M0) and 
likely formed most of their stellar mass at z > 5. To a limit of 0.5 x 10 M their number density is ~ lOx 
lower than that of local early-type galaxies. Furthermore, we use the IRAC photometry to derive rest-frame 
near-infrared J,H, and K fluxes. The DRGs and LBGs together show a large variation (a factor of 6) in the 
rest-frame ^-band mass-to-light ratios (M/Lk), implying that even a Spitzer 8^m-selected sample would be 
very different from a mass-selected sample. The average M/Lk of the DRGs is about three times higher than 
that of the LBGs, and DRGs dominate the high-mass end. The M/Lk ratios and ages of the two samples appear 
to correlate with derived stellar mass, with the most massive galaxies being the oldest and having the highest 
mass-to-light ratios, similar as found in the low-redshift universe. 
Subject headings: galaxies: evolution — galaxies: high redshift — infrared: galaxies 



1. INTRODUCTION 

Galaxies at z > 2 exhibit very diverse properties: they range 
from the blue Lyman Break Galaxies (LBGs) which are bright 
in the rest-frame ultraviolet (Steidel et al 1996a,b) to the Dis- 
tant Red Galaxies (DRGs) which are generally faint in the 
rest-frame UV and have fairly red rest-frame optical colors. 
The DRGs are selected in the observers' near-infrared (NIR) 
by the simple criterion J s -K s > 2.3 (Franx et al. 2003, van 
Dokkum et al. 2003). The variety in the galaxy population at 
z > 2 is comparable to that seen in the local universe, where 
colors range from very blue for young starbursting galaxies 
to very red for old elliptical galaxies. An urgent question is 
what causes the red colors of DRGs at z > 2. Are they "old 
and dead", or are they actively star forming and more dusty 
than U-dropout galaxies? 

First analyses of the optical-to-NIR Spectral Energy Dis- 
tributions (SEDs) and spectra have suggested that both ef- 
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fects play a role: they have higher ages, contain more dust, 
and have higher mass-to-light (M/L) ratios in the rest-frame 
optical than LBGs (Franx et al. 2003, van Dokkum et al. 
2004, Forster Schreiber et al. 2004; F04). Furthermore, many 
have high star formation rates > 100 M Q yr _1 (van Dokkum 
et al. 2004, F04). Unfortunately, the number of DRGs with 
rest-frame optical spectroscopy is very small. Inferences in- 
evitably depend on S EP fittin g, which has large u ncertainties 
(see e.g., Papovich, Dickinson, & Ferguson 2001). The SED 
constraints on the stellar and dust content are expected to im- 
prove significantly by extending the photometry to the rest- 
frame near-infrared. 

Here we present the first results on rest-frame NIR pho- 
tometry of DRGs and LBGs in the Hubble Deep Field South 
(HDFS) as observed with the Infrared Array Camera (IRAC; 
Fazio et al. 2004) on the Spitzer Space Telescope. Where 
necessary, we assume an Hm = 0.3, J7a = 0.7, cosmology with 
Ho = 70 km s _1 Mpc _1 , and a Salpeter Initial Mass Function 
(IMF) between 0. 1 and 100 M Q . Magnitudes are expressed in 
the AB photometric system unless explicitly stated otherwise. 

2. THE OBSERVATIONS, PHOTOMETRY, AND SAMPLE SELECTION 

We have observed the 5 arcmin 2 HDFS/WFPC2 field with 
the IRAC camera integrating 1 hour each in the mid-infrared 
(MIR) 3.6, 4.5, 5.8, and 8/im-bands. The observations, data 
reduction, and photometry will be described by Labbe et al. 
(in preparation). Briefly, we reduc ed and calibrated th e data 
using standard procedures (e.g., Barmbv et al. 2004). The 
limiting depths at 3.6, 4.5, 5.8, and 8/zm are 25.0, 24.6, 22.6, 
and 22.5 mag, respectively (5rr, 3" diameter aperture). We 
PSF-matched the IRAC images and a very deep Zf s -image 
(Labbe et al. 2003) to the 8/im-band. Nearby sources were 
fitted and subtracted to avoid confusion (see Labbe et al. 
2004). We measured fluxes in 4."4 diameter apertures. The 
K S -IRAC colors were then combined with the HDFS-catalog 
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of Labbe et al. (2003) resulting in 1 1 -band photometry from 
0.3 to 8/im. In addition to the photometric errors, we add in 
quadrature an error of 10% to reflect absolute calibration un- 
certainties of IRAQ 

We selected DRGs in the field of the HDFS using the crite- 
ria of lFranx et alJ 1120 03) yielding 14 galaxies with photomet- 
ric redshifts (Rudnick et al. 2003) in the range 1 .9 < z < 3.8. 
One blended source was excluded. In addition, we selected 
isolated LBGs in the same field to the same ^-band limit 
and in the same red shift range, using the U-dr opout criteria of 
i Madau et alJ {l996) on the WFPC2 imaging dCasertano et alJ 
2000). The two samples are complementary, with only 1 
source in common. For both samples, the mean redshift is 
z = 2.6 with a standard deviation of 0.5 11 . 




-0.5 0.0 0.5 1.0 1.5 2.0 
K s - 4.5/i.m 

FIG. 1. — Observed I — K s versus K s — 4.5/im color-color diagram of two 
samples of z > 2 galaxies. Distant Red Galaxies (DRGs; filled circles and 
stars), occupy a different color region than z, ~ 2.5 Lyman Break Galaxies 
(LBGs; diamonds). DRGs that are best described by an old single-age burst 
model (SSP; stars) have colors distinct from those better described by con- 
stant star-forming models (CSF) and dust, indicating that IRAC fluxes may 
help in separating these populations. The curves show the color-evolution 
tracks of Bruzual & Chariot (2003) models at a fixed z = 2.6: an SSP model 
with ages ranging from 0.3 to 3 Gyr (dashed line) and two CSF models with 
ages ranging from 0.1 to 3 Gyr, and reddenings of Ay = 0.6 and Ay = 2.0 
respectively (solid lines). The vector indicates a reddening of Ay = 1 mag for 
a Calzetti et al. (2000) law. 

3. RED GALAXIES AT Z > 2: OLD OR DUSTY OR BOTH? 

In Fig. [2 we analyze the I — K s versus the K s -4.5fj,m col- 
ors. The K-4-.5fj.rn color has a sufficiently large baseline 
and much higher signal-to-noise ratio than the K s -5.8fim or 
K s -8/im colors. The DRGs and LBGs lie in very differ- 
ent regions, with little overlap. The mean K-4-.5fj.rn color 
of the DRGs is significantly redder than that of LBGs, con- 
firming that DRGs have higher M/L ratios in the rest-frame 
optical. The red I — K„ and K-4.5fim colors of the DRGs 
imply that they must be prominent in IRAC selected sam- 
ples. Indeed, the majority of red z-3.6/im sources selected 
by Yan et al. (2004) in the GOODS survey satisfy the DRG 
selection criteria. We show color tracks for stellar popula- 

11 The photometric redshift accuracy of DRGs from other studies is 
\zspec~ z p hot\/(l +z S pec) ~ 0.1 (F04, Wuyts et al. in preparation). The esti- 
mate is based on direct comparison of 16 galaxies with both photometric and 
spectroscopic redshifts. Uncertainties in rest-frame magnitudes and model 
parameters are based on Monte Carlo simulations that take into account un- 
certainties in the observed fluxes and redshift estimates. 



tion models (Bruzual & Chariot 2003; BC03), redshifted to 
a fixed z = 2.6: a single-age burst model (SSP) and two mod- 
els with constant star formation (CSF), and reddenings typical 
for LBGs (A v = 0.6, Shapley et al 2001), and typical for DRGs 
(Ay = 2.0, F04). The K s -4.5fim colors of LBGs are consis- 
tent with low reddening models, whereas most of the DRGs 
lie in the area of models with substantial extinction. 

3.1. Dead galaxies at z = 2 — 3 

Three of 13 DRGs (indicated by star symbols in Fig. 1) lie 
well outside the area of constant star-forming models. Their 
I—K s colors are too red for their K— 4.5 /mi color and they lie 
close to the line of an old single-age burst model with ages 
of 2 - 3 Gyr. The candidate old, passive galaxies are shown 
in Fig. |2 Their SEDs are very well fit by an old single burst 
population. The gray model curves are predictions based on 
fits to the exceptionally deep optical/NIR data only, demon- 
strating that the Spitzer fluxes lie very close to the expected 
values. Dead galaxies have blue 4.5/im colors and with 
Spitzer we can now effectively separate them from dusty star- 
forming DRGs, which are red in K-4.5fj.rn. 

Models with ongoing star-formation and dust-reddening fit 
the SEDs badly. Star-formation histories (SFHs) with expo- 
nentially declining star-formation rates (SFRs) and timescales 
> 300 Myr can be ruled out at the 99.9% confidence level. 
For the marginally acceptable models, the ratio of on-going 
to past-average SFR is SFR(t)/ < SFR >= 0.001 indicating 
that these 3 galaxies are truly "red and dead". The best-fit 
ages and implied formation redshifts depend on the assumed 
metallicity Z, as expected. Super-solar Z = 0.05 models give 
a mean age of < t >= 1 .4 Gyr and mean formation redshift 

< Zf >= 5, solar Z = 0.02 models yield < t >= 2.9 Gyr and 

< Zf >S> 5, while sub-solar metallicity models fail to provide 
good fits to the IRAC fluxes as they are too blue. Hence we in- 
fer from the models that the "dead" galaxies had formed most 
of their stellar mass by z ~ 5 (Z = 0.05) or much earlier z S> 5 
(Z = 0.02). 

The number density of "dead" galaxies is 1.9 
xl0~ 4 /i7 Mpc" 3 , assuming a top hat redshift distribu- 
tion between z = 2 and 3. For stellar masses > 0.5 x 10 n M Q 
and the same IMF, that number density is 10 x lower than 
that of early-type galaxies in the nearby universe (Bell et 
al. 2003). Obviously, the sampled volume is too small to 
allow firm conclusions, but the result is indicative of strong 
evolution of passive galaxies from z = 2.5 to z = 0. 

3.2. Dusty star forming galaxies 

In Fig. |3 we show the SEDs for 3 of the 8 DRGs whose 
optical/NIR fluxes were better fit with CSF models and red- 
dening by a Calzetti et al. (2000) dust law. The observed 
MIR flux points are often different from predictions based on 
shorter wavelengths, especially where the flux density F\ was 
still rising at the observed K-b&nd. Hence, for the very dusty 
DRGs (Ay > 1 .5) the MIR fluxes help to better constrain the 
age and dust in the models. The average ages and extinction 
of the fits changed mildly after inclusion of the IRAC data: 
from < t >=1 .1 to 1.3 Gyr, and <A V >=1 .9 to 1.5 mag, re- 
spectively. 

Finally, two remaining galaxies could not be fit well by 
any model, most likely due to emission line contamination. 
One of them is spectroscopically confirmed as a very strong 
Lyman-a emitter (Wuyts et al. in preparation). However, the 
generally good fits to the UV-to-8/im SEDs of the DRGs in- 
dicate that the red J s - K s colors were caused by old age and 
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FIG. 2. — The UV-to-8/im spectral energy distributions (SEDs) of Distant Red Galaxies (DRGs). The best-fit Bruzual & Chariot (2003) models to the full 
SEDs (black solid lines) are shown, but also predictions based on the ultradeep optical and near-infrared fluxes only (gray solid lines). Top row: the 3 DRGs 
best fit with an single-age burst (SSP) model. The mid-infrared (MIR) IRAC fluxes (filled circles) directly confirm the predictions. Old and dead DRGs have 
bluer K s — 4.5/mi colors than very dusty star-forming DRGs. Galaxy 767 shows a flux excess at 8/rni, possibly related to AGN activity. Bottom row: the MIR 
predictions for galaxies better fit with constant star forming models and Calzetti et al. (2000) dust-reddening. For highly-reddened galaxies, the MIR observations 
can be very different from the predictions. Here IRAC fluxes help to better constrain the total age, dust content, and stellar mass in the models. 



dust-reddened star light, and not other anomalies. We note 
that one of the dead galaxies (see Fig. 1) has an apparent flux 
excess at 8^m, suggesting the presence of an obscured AGN 
which st arts to dominate the flux in the rest-frame /T-band 
(see e.g.. lSternetalJ2005l) . 

4. MASS-TO-LIGHT RATIO VARIATIONS AT Z = 2 - 3 

We next study the M/L ratios inferred from the SED fits 
described above, adopting the best fitting SFH (SSP or CSF). 
Figure 3 a shows the modeled M/Lk (rest-frame K) versus 
rest-frame U — V . The curves show the dust-free BC03 model 
M/Lk ratios. Galaxies generally lie to the red of the model 
curves as a result of dust attenuation. The fits imply a large 
range in M /Lk for DRGs and LBGs together (a factor of 6). 
Furthermore, the M/Lk ratio correlates with rest-frame U — V 
color, where galaxies red in the rest-frame optical have higher 
M/Lk ratios than blue galaxies. The average M/Lk of the 
DRGs is 0.33 ±0.04, about three times higher than that of the 
LBGs. Figure 3b shows the derived M/Lk against the derived 
age. The M/Lk correlates well with the age, as expected from 
the models. The role of extinction is greatly reduced in the 
rest-frame TT-band, implying that the differences in the M/Lk 
ratios for DRGs and LBGs are driven by age differences. 

Figure 3c shows the M/Lk ratio against stellar mass M*. 
DRGs dominate the high-mass end. The highest-mass galax- 
ies (> 0.5 x 10 n M Q ) in this sample all have high M/L K ra- 
tios, and here the M/Lk ratio does not depend strongly on M* . 
At lower masses galaxies have much lower M/Lk ratios. This 
may be a selection effect caused by our magnitude cut-off, 
as we would miss low-mass galaxies with high M/Lk- How- 
ever, high-mass galaxies with low M/Lk would be detected 
if they existed, hence the lower envelope of the M/Lk versus 
M* distribution is real. Two intriguing possibilities are that 
the mean M/Lk decreases towards lower stellar mass, or that 
the intrinsic scatter in M/Lk increases. 

Our DRGs and LBGs were selected in the observed 
A^-band (rest-frame V-band). We find a large variation in 



M/Ly (a factor of 25), hence selection effects play a major 
role in NIR studies at high redshift. Selection using MIR ob- 
servations with Spitzer would improve the situation, but the 
wide range in M/Lk ratios and ages found here indicates that 
even a IRAC 8/im-selected sample would still be very differ- 
ent from a mass-selected sample. 

5. DISCUSSION AND CONCLUSIONS 

We have presented rest-frame UV-to-NIR photometry of a 
sample of DRGs and LBGs at z = 2- 3 in the HDFS. These 
galaxies span a wide range in properties, similar to low- 
redshift galaxies. The rest-frame NIR photometry from IRAC 
helps significantly: first, by allowing us to separate "old and 
dead" from dusty star forming DRGs using only the observed 
I,K S , and 4.5/im-band, and second, by improving model con- 
straints on the heavily obscured DRGs. 

The wide range in galaxy properties at z = 2.5 raises sev- 
eral important issues. First, it demonstrates it is impossible 
to obtain mass-selected samples photometrically. Even in the 
rest-frame /T-band, the M/L ratio varies by a factor of 6 for 
the DRGs and LBGs in the HDFS. Second, we need to under- 
stand what produces these variations. If a relation between to- 
tal stellar mass and M/Lk exists (see Fig. 3c) then stellar mass 
might be driving the variations. Deeper IRAC data is needed 
to establish this well, as incompleteness effects may play a 
role. It is tantalizing that Kauffmann et al. (2003) find a simi- 
lar correlation at z = 0. 1 , with the most massive galaxies being 
the oldest and having the highest M/L ratios. The authors also 
found that above a stellar mass M* = 6 x 10 10 M Q galaxy prop- 
erties correlate only weakly with M„, similar to what we find 
at z = 2 - 3. The simplest explanation is that we observe the 
same galaxies at z = 0. 1 and z = 2.5, although we note that hi- 
erarchical formation scenarios predict significant merging for 
galaxies at z > 2. Alternatively, we observe merging processes 
occurring at a critical mass of about 6 x 10 I0 M Q . A partial ex- 
planation may be a simple relation between mass and metal- 
licity: higher mass galaxies might have higher metallicity, and 
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FIG. 3. — The mass-to-light ratios in the rest-frame S'-band (M/Lk) of distant red galaxies (DRGs) and z = 2 — 3 Lyman break galaxies (LBGs). Symbols 
are as in Fig. 1. (a) M/Lk versus rest-frame U — V color, showing a clear correlation. On top we indicate the corresponding observed J s — K s color and the 
Js-Ks > 2.3 limit (upward arrow). Color tracks of an SSP (dotted line) and CSF model (dashed line) are shown, (b) The relation between M/Lk and best-fit 
age. The M/Lk is more sensitive to the age of the stellar population than to dust extinction, hence the DRGs are on average older, not just more dusty versions of 
LBGs. (c) M/Lk versus stellar mass. The highest-mass galaxies all have high M/Lk ratios. The mass completeness limit (downward arrow) corresponds to our 
observed K SJOt = 22.5 magnitude cut-off and a maximal M/L. A selection by rest-frame K-light corresponds to taking everything to the right of the constant Lk 
line (dashed line). Obviously this is still far from a vertical cut-off which would represent a selection by stellar mass. 



thus be redder. 

It is very noteworthy that there are "red and dead" galaxies 
at z = 2-3. Previous authors found such galax ies at z = 1-2 
(e.g.. lCimatti et alJ2004tlMcCa rthv et a l 120041) . Whereas the 
number density of "dead" galaxies at z = 2.5 is probably much 
lower than at z = 0, the mere existence of these systems at such 
high redshift raises the question what caused such a rapid and 
early decline in star formation. Our model fits imply that their 
star formation stopped by z = 5 or higher, close to or during the 
epoch of reionization. Possibly a strong feedback mechanism 
caused this, such as an AGN or galactic-scale outflow. We 
note that these galaxies are among the most massive galaxies 
at z = 2.5, and hence were probably at the extreme tail of the 
mass function at z = 5. 



Obviously, many questions remain unanswered by this 
study. Only a very small field has been studied, and simi- 
lar studies on wider fields are necessary. Finally, the M/L 
ratios derived here remain model dependent and vary with the 
assumed SFH, metallicity, and IMF (e.g., Bell et al. 2003). 
Direct mass determinations are presently very challenging, 
but are becoming more important to understand the z = 2 - 3 
galaxy population. 
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